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SECTION 1 


BACKGROUND 


In 1970, AFGL contracted DBA Systems to perform a study 
of various reduction approaches to satellite altimetric measurements 
of the oceanic geoid surface. This study was initiated in view of 
the then upcoming launch of the GEOS-3 Satellite which was designed 
to carry a radar altimeter. The primary objective was to derive 
alternative reduction approaches that may remove the necessity for 
extremely accurate reference orbits of GEOS-3 required by convent- 
ional approaches. 

• i 

Conventional approaches to the reduction of satellite 
altimeter data utilize long arc orbit integration that requires 
extensive tracking from ground based trackers to maintain extremely 
accurate orbits, A widely held premise was that in order to exploit the 
satellite altimeter of 1 meter accuracy for geoid improvement, it was 
necessary that the radial component of satellite position be known 
to better then + 1 meter. An example is extracted from reference 11 
that states "Because the orbit will be used in combination with 
altimeter measurements in the data reduction process leading to geoid 
improvement, it is necessary that the accuracy in S/C height as cal- 
culated from the orbit be known to better than + 1 meter". Problems 
arose that made the above requirement very difficult. Among these 
problems are (a) unresolved biases in a given tracker could induce 



looat'ized systematic errors in the confuted orbit which could he 
transferred to the local geoid^ (b) vnmodeled "perturbations caused 
by drag radiation pressure, (c) errors in long arc orbits result 
from integrated effects of errors in the adopted geopotential function 
and (d) most important is the extensive effort and cost of global 
tracking networks and computer processing required to provide such 
accurate reference orbits. 

In view of the above difficulties, our investigation was 
directed toward alternative approaches that would fully exploit the 
1 meter altimeter accuracy for geoid determination with less stringent 
requirements for orbital accuracies. This led to the investigation 
of the feasibility of utilizing the short arc technology that AF6L/ 

DBA had developed in previous satellite geodesy programs. Such an 
approach was envisioned to involve the simultaneous recovery of 
state vectors defining as many as several thousand short arcs (sub- 
ject to weak or prior constraints) along with the mathematical model 
used to define the geoid surface. The estimated orbit accuracy 
requirement was + 100 meters, which is easily attainable from 
routine global tracking. 

The short arc approach is defined as orbital arcs no 
greater than one fourth revolution. It was shown in Brown {1967) 
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that accuracies of integration of better than 1 meter can be attained 
in short arc reductions provided that 

a) sphen^cdl hcmnoni-es at least through 
(n^m) = (4,4) are exeroi-sed in the 
integration 

and 

b) all six orbital -parameters 
(Xo i Yo » Zo 5 Xo» Yo 5 Zo) 

at mid-aro are free to adjust to 
best aoeormodate the aotuat orbit. 
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SECTION 2 


INTRODUCTION OF THE SHORT ARC APPLICATION 
TO REGIONAL GEOID REDUCTIONS 

1 

The initial investigation established that the short arc 
approach was feasible and offered a flexible means of utilizing 
altimeter data for detennining the geoid surface for a wide range 
of applications. The initial consideration of the short arc method 
employed spherical harmonics to mathematically define the geoid 
surface. This development is well suited to global representations 
of the geoid undulations and has been implemented into a computer 
program {SAGG) for the combination of Satellite Altimetry and 
Ground Gravity reduction. However, it was anticipated that the 
early data collection phase of the GE0S~3 program would concentrate 
on altimeter measurements over the North Atlantic calibration area. 
Additionally, it was felt that investigations of fine detail over 
limited localized regions would require an extravagant spherical 
harmonic expression. Therefore, alternative means for analytical 
representation of the geoid surface were explored. The most 
suitable model for the proposed application was a derivation of the 
spheroidal multiquadric analysis developed by Ear‘c3y (1978). In 
Hardy's exercise of the multiquadric analysis, the nodes (aoeffieient 
Gomputation point) correspond to observed data points. In our deriva- 
tion, the nodes do not necessarily correspond to data points, are 
relatively limited in number and, at the outset, are evenly distributed. 
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The primary objective of this project was the short arc 
determination of the North Atlantic geoid, utilizing the spheroidal 
multiquadric analysis. This led to the development of the computer 
program SARRA {Short Are Eeduation of Radocr Attimetry) , It was 
envisioned that the altimeter observation would cover the entire 
North Atlantic Ocean at spacings of approximately a 1° by 1“ grid 
with an initial selection of nodal points at 5° by 5® grid as 
illustrated by Figure 1. 

As data became available, it was evident that there was 
a high concentration of the observations over the calibration area with 
sparse measurements over the northeastern portion of the North 
Atlantic Ocean and almost none in the southeastern portion. The 
results of the reduction indicated a close detailed agreement of 
the dense data area with other geoid models {Marshy Vincent ^ Strange) ^ 
but the geoid in the. sparse areas {northeastern) seemed tilted toward 
the more dense areas. This geoid behavior caused concern with the use of 
the multi quadric analysis when the data is non-uniform and sparse in cer- 
tain areas. Consequently, other possible surface models were investigated. 

Upon suggestions by Dr. Donald Eckhart (AFGD* and extensive 
theoretical investigation by Dr. George Blaha (DBA)**^ the use of the 
covariance function was introduced as a second option to the computer 
program SARRA. The two options now offer flexibility depending on 
the characteristics of the available data. Both options are demonstrated 

* Air Force Geophysics Laboratory, LGBanscom Air Force Base, Ma. 

Personal communication. 

** DBA Systems, Inc., Melbourne, FL. Personal communication. 
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.FIGURE 1. Illustrating ground tracks of every fifth pass of 
approximately 200 passes of GEOS C considered in 
' hypothetical short arc determination of fine 
structure of North Atlantic geoid. Indicated 
tracks generate approximate cells. Selected 

nodes for initial representation of North Atlantic 
geoid by means of spheroidal multi quadric functions 
are located nominally at alternate corners of^ 

5“x5° cells as indicated by solid dots. Initial 
set of approximately 150 nodes is subsequently to - 
be augmented by additional nodes at locations 
indicated by residuals from initial reductions. 
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in the results portion of this report. Details of the multiquadric 
analysis and the covariance function models are presented in Section 
4.1 of this report. The covariance function that is implemented into 
SARRA refers to a chosen reference ellipsoid. 

The advantage of the covariance function is realized 
when data is irregularly distributed with some surface areas 
being very sparse in measurements. The covariance function 
implies some apriori knowledge of the geoid behavior through 
the spherical harmonic coefficients. The spheroidal multiquadric 
analysis is completely dependent on observation data and reproduces 
a surface model as the best least squares fit of the observation 
data. The advantage of using the multiquadric analysis model is 
realized when the area to be processed is covered with an adequately 
dense set of altimeter measurements. One unique feature of the multi- 
quadric analysis is that there is essentially no limit on the size 
of the area to be processed and, provided enough data, one can 
obtain as much detail as the data provides by the proper selection 
of nodal points. 
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SECTION 3 


INVESTIGATIVE EXPERIMENTS CONDUCTED TO DETERMINE 
THE FEASIBILITY OF THE SHORT ARC APPROACH 


In the development of the computer program SARRA, a series 
of investigations were conducted for the purposes of determining 

a) the has'lo feasih-iUty of the appZiaation of 
the shovt ca?e teahnotogy to geoid recovery 
from sateZZtte aZtimetrys 

b) the most suitahZe math models for represent— 
i-ng the geoZd surface Zn regard to tZie Zn- 
tended ccppZZoatZony 

c) testZng the overaZZ concept wZth sZmulatZon 
for Zts uZtZmate accuracy potentdaZ, 

d) formuZatZon of the overaZZ concept Znto a 
reduction program capabZe of routZneZy 
processZng the observatZon data when 
GEOS-3 became operational 

and 

e) deveZop auxZZZary programs for the pre and 
post processing required in an operatZonaZ 
environment. 

The key steps in this investigation are summarized in this section. 
However, the full detail may be found in references 1 and 11. 

It was surmised in Brown {l97d) that the spheroidal 
multi quadratic model would be well suited to detailed representa- 
tion of the geoidal surface over such local regions as the North 
Atlantic {where intensive testing of GEOS-C was to be conducted). 
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In such applications, the model has the virtue of relative simplicity 
and is especially attractive in the flexibility afforded by the pro- 
cess referred to as nodal dens-tfioation. This is an adaptive process 
wherein the original set of regularly spaced nodes {defining the 
multiquadrio fimotion) are supplemented locally as needed to improve 
the fit over irregular areas inadequately modeled at the outset. 

Despite its theoretical attractiveness, the spheroidal 
multiquadric model had to be implemented and tested. The next section 
outlines the measures taken in evaluating the model under a well 
controlled experiment with data typical of geoid features. 

3.1 Preliminary Exercises of the Spheroidal Multi quadric Model 

In order to ascertain the basic adequacy of the spheroidal- 
multi quadric model, two preliminary numerical tests were performed. 
The first test was to determine if the model was inherently capable 
of providing a reasonably good fit to the geoid; the second was to 
test the model as incorporated in SARRA with extensive simulations. 
These simulations will be presented in the next section(Section 3.2). 

The first test used as a data base, the gravimetric geoid 
of a portion of the North Atlantic produced by Marshy Strange and 
Vincent {l9?2). A sample of about 250 spot elevations of the 
geoid were extracted from the contour map of this geoid. These 
geoidal heights were treated as if they were direct observations. 
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The multiquadric model was then fitted by least squares to this data 
set using an initial set of 112 nodes spaced at five degree intervals 
(see Figure 2). The process of nodal densifi cation was invoked in 
three iterations of the adjustment to introduce fresh nodes in re- 
gions leading initially to residuals in excess of two meters. This 
led ultimately to the incorporation of 22 additional nodes, several 
of which were in the region of the severe undulation over the Puerto 
Rico Trench. The final rms error of the fit of the multi quadric 
model turned out to be an altogether acceptable 0.75 m py:> 6 < 2 ision 
of the dig-itfzation of the oontour map was deemed to he not much 
better than 0.5 m) . This test established two things: 

(a) that the spheroidal multiquadratie function 
could provide an aocunate representation of 
regional undulations of the geoid^ 

and 

(b) that the process of nodal densification 
provides a validj effective and efficient 
means for extending the model as needed 
for local improvement of fit. 

Figure 3 is the contour produced by the reduction and Figure 4 
represents contour of closeness of fit. 

With the fundamental soundness of the spheroidal multi - 
quadric model firmly established as a result of the first preliminary 
test, steps were taken to incorporate the model into SARRA. After 
the revised program had been checked out on DBA's Xerox Sigma 5 
computer, a small scale simulation was performed to provide a 
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80’ 75’ 60* 45’ 30’ 15* 

• Supplemcniar/ nodes 
-t- Supplementary Data 

The Initial node selection is represented by the intersection of the grid lines of figure. 

_ The initiol doto wos interpolated at each 2* grid starting at latitude 17 ond longitude 15 west 
through loti,tude 44’ and longitude 79 * west. 


FIGURE 2. Map Showing Contours of Gravimetric Geoid and Locations 
of Selected Nodes Employed in Preliminary Testing of 
the Spheroidal Multi quadric Model 





FIGURE 3. RMS Fit = .75 meters 

Grid Level = 2 degrees 

Nodal Density = 10 degrees 








tentative indication of accuracies to be expected in applying the model to 
the determination of geoidal undulations over the North Atlantic. 

For this purpose, a set of 40 passes was generated, half of which 
were ascending (SE-NW) and half of which were descending (NE-SW). 

Nodes (a total of 7l) were placed at alternate points on a five degree 
grid and a low sampling rate of one point ever 70 seconds was adopted 
[this corresponds to a spacing of about five degrees along each arc), 

A priori one sigma constraints of 5m were exercised on the state vectors 
for each of the 40 arcs and a standard deviation of Im was assigned to 
the altimetry. Under these assumptions, the error propagation performed 
by SARRA indicated that rms errors in the recovered geoid could be 
expected to range from a low of 0.67 in the interior of the region, 
to a high of 1.93m near the boundary. This result is not considered 
to be particularly realistic because of the obviously inadequate num- ' 
ber of nodes exercised in the model. It does, on the other hand, 
provide some indication of the potential power of the method, parti- 
cularly when consideration is given to the very low sampling rate 
that was adopted. 

3.2 Computer Simulation of the Recovery of the Geoid Undulation 

Over the North Atlantic 

The successful outcome of the preliminary test made it logical 
to proceed with large scale simulations that would be more nearly repre- 
sentative of what could be expected from the GEOS-3 tests to be conducted 
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over the North Atlantic. For this purpose, the revised version of 
SARRA was converted to run on the CDC-6600 computer at AFGL. The 
simulations performed, consisted of two runs based on sets of 320 
passes. In all simulations, the nodes were spaced at regular five 
degree intervals, the total number being 126. In the case of the 
200-pass reductions, the altimeter observations were assumed to be 
spaced at nominal one degree intervals along each arc [tMs aorres- 
ponds to a time interval of 14 seconds). In the case of the 320- 
pass reduction, the chosen spacing corresponds to 0.625 degrees, or 
a time interval of 8.5 seconds. Schematic layouts of the ground 
tracks for the two sets of simulations are presented in Figures 5 
and 6. Other pertinent data concerning the distinguishing assump- 
tions underlying the various simulations are indicated in Table 1. 

In Cases 3, 4, 5 and 6, the sigma adopted for altimetry, 
namely 0,37m,. is representative of what would be expected from a 
data compaction process based on the following considerations: 

(a) an original sampling, rate of tm per second, 

(b) an original sigma of 1.0m, and 

(c) use of a third order midpoint, cubic 
polynomial filter employing 17 points. 

In Cases 1 .and 2, no data compaction scheme was considered to have 

been exercised; accordingly, here the expected sigma of the raw 

altimeter observations was exercised at the aforementioned data rate 

of one point every 14 seconds. 
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The output of primary interest from thp simulations consists 
of the- propagation of the specified errors in the input into the com- 
puted geoidal heights. Such error propagations were, performed for 
geoidal points corresponding to the nadirs of the data points. In 
due course, such results will be employed to generate contour maps 
defining the variation. in the error of the recovered geoid. In this 
section, consideration is limited to the results indicated in 
Table 2; here, the high, low and average values of the standard 
deviations of the recovered geoid are listed for each of the Simula^ 
ti ons . 


TABLE 1. Basic Assumptions Underlying Various Simulations 


■■ 

NUMBER 
OF PASSES 

A PRIORI ONE SIGMA 
ORBITAL CONSTRAINTS 
(Position/Velocity) 

ALTIMETRYcj(m} 

1 

200 

20m/. 02 m sec ^ 

1.0 

2 

200 

15m/. 015 m sec ^ 

1.0 

3 

320 

■15m/. 015 m sec"^ 

0.37 

4 

320 

-1 

15m/. 015 m sec 

0.37 

5 

320 

5m/. 005 m/sec ^ 

0.37 

6 

320 

-1 

Im/.OOl m/sec 

0.37 
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Table 2. Key Results of Simulations 


Case 

Number 
of Passes 

One Sigma 

Orbital Constraints 
(Posi ti on/Vel oci ty ) 

Sigmas of Recovered Geoid 


Low 

mmm 

1 

200 

20m/. 02 m sec"^ 

1.52 

1.34 

1,42 

■ 

200 

15m/. 015 m sec"^ 

.96 

.82 

.89 

B 

320 

20m/. 02 m sec“^ 

1.03 

.97 

1.00 

B 

320 

15m/. 01 5 m sec”^ 

.54 

.46 

.51 

B 

320 

5m/. 005 m sec"l 

.33 

.22 

.27 

H 

320 

Im/.OOlm sec"^ 

.13 

.06 

.09 


The results indicate that general accuracies comfortably 
better than Im {rms) are to be expected from the reduction of 320 
passes, subject to a priori orbital constraints on the order of 
15m (Case 4). Because orbital accuracies of considerably better 
than 5m are a reasonable expectation for the precise reference orbits 
ultimately to be generated, the results from Cases 5 and 6 suggest that 
geoidal accuracies on the order of a few tenths of a meter are potent- 
ially attainable through the application of the Sho:ct Arc Method. 
Attainment of such accuracies in practice will, of course, entail the 
appropriate application of nodal densifi cation. This, in turn, will 
require processing of additional passes to maintain a specified level 
of accuracy. A reasonably conservative extrapolation, in our view, 
is that with the exercise of some 50 to 75 well-placed, additional 
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nodes (raising the totaU from 275 to 200) and with the incorporation 
of perhaps an extra 100 passes [for a total of onev 400) , one can 
expect SARRA to produce geoidal accuracies generally better than 0.5m 
when precise reference orbits are employed in the adjustment with one 
sigma a priori constraints of 5m {or better). 

3.3 Reduction of the First Available Satellite Altimetry 

Data from Skylab 

The first satellite altimetric data for experimental process- 
ing with the SARRA computer program became available from altimetric 
observations made by Skylab [SL~2) . Two passes were processed through 
the SARRA computer program for the purpose of testing the basic data 
flow procedures and to evaluate the adequacy of the program in reducing 
the data. Figure 7 shows the ground track of the two passes plotted by 
the DBA CALCOMP plotter. Pass 9 begins off the east coast of the United 
States and extends approximately 15 to 20 degrees below the equator. 

Pass 4 begins off the east coast and extends over Puerto Rico. 

Mid-arc state vectors were computed from ephemeris data 
provided by NASA for initial orbital estimates of each arc. Data 
editing procedures were designed and performed on the observation 
data prior to SARRA processing. These included, primarily, an 
automatic rejection of gross error indicated by the residual altimetry 
parameter in the SL-2 EREP format and a three-sigma criterion edit 
after a ninth order polynomial data fit. 
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FIGURE 7. Skylab Ground Tracks 


The approach taken In this experiment was to process the 
passes through SARRA and examine the behavior of the minimized resid- 
uals. Since Pass Four extended over the Puerto Rico Trench vicinity, 
the residuals were examined particularly for undulations in that 
area. The residuals after SARRA adjustment are shown in Figure 8. 

A sudden drop in the residuals is immediately noticed 126 seconds 
from epoch. This section of the pass was extracted and processed 
again with only the observation occurring over this region. The 
residuals after adjustment are shown in Figure 9, This residual 
profile was plotted on a figure {fi.gio'e lO) reproduced from previous 
comparison of SL-2 versus gravimetric geoid [Vincent, Strange, Marsh). 

The residuals of Pass Nine are presented in Figure 11. 

The residuals show no abrupt changes that would indicate any possible 
sharp geoid undulation. There is a gradual but smooth variation 
throughout the pass of approximately 7.0 meters that may be due to a 
combination of geoid undulation and sea state variations. 

The importance of this experiment was simply to gain exper- 
ience with the SARRA program with real observation data and to illus- 
trate the sensitivity of the program in- detecting geoid undulations. 
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Residuals (Meiers) 



31.0 Lot. 25.5 Lot. 24.5 Lot. 19.9 Lot. ia.4Lot. 

8C.2 Lon. 74. 1 Lon. 73.0 Lon. 48.7 Lon. 47.4Lon. 


Latitude and West Longitude 


FIGURE 8. Skylab Altimeter Data Pass #4 



Residuals (Meters) 



FIGURE 9. Skylab Altimeter Data Pass #4 
Puerto Rico Trench Only 
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Result from DBA SKYLAB SL-2 altimeter data plotted on a reproduction of Figure 7 from Marsh. Vincent 
"DETAILED GRAVIMETRIC GEOID FOR THE GEOS-C ALTIMETER CALIBRATION AREA". 


FIGURE 10. Comparison of Gravimetric Geoid with Skylab Altimeter Data 





3.4 


Preliminar.y Reduction of GEOS-3 Altimetr.y Data Over 
the North Atlantic Ocean 


Some preliminary reductions were performed with GEOS-3 
altimetric observations concentrated over the western portion of 
the North Atlantic Ocean. The observational data was provided to 
AFGL by NASA and consisted of 112 passes {Figure 12). These data 
sets provided an opportunity to test the basic mathematical models 
used to represent the geoid surface. The results of these reduc- 
tions were presented at the American Geophysical Union {agu) fall 
meeting held in December 1976 at San Francisco. 

The distribution of this data set was not adequate for a 
geoid undulation determination of the entire North Atlantic and was 
not suitably located for the detailed examination of such fine fea- 
tures as the Puerto Rico Trench. However, it did provide adequate 
data to evaluate the performance of SARRA computer program in an 
operational environment and demonstrate the potential accuracy of 
the approach. 

The short arc approach to these reductions may be visualized 
along the following lines: 

(a) reference orbits accurate to approximately 
20 meters may he obtained from routine 
global tracking; 

(b) the reference orbits are divided into sub- 
arcs ^ situated over oceanic regions and are 
limited in length to no more than a quarter 
revolution; 
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(c) eaoh sub-aro is treated as an independent orbit 
with the epoah at mid-are having a state vector 
subject to a priori constraints consistent with 
the estimated accuracy of the reference orbits; 

(d) observational equations are formed from satellite 
altimetrie measurements for each sub-arc for the 
mathematical model chosen to represent the oceanic 
geoid; 

(e) the adjustment simultaneously recovers the coef- 
ficients of the geoid surface model and revised 
estimates of the state vectors for all sub-arcs 
(whlah 0A.2. anlbmittd -in numbeAj. 


3.4.1 Unknowns and Constraints Used in the Reductions. The SARRA 
program is designed to determine surface coefficients as designated 
by nodal points that lie within the boundaries of a dense selection 
of height measurements. In addition to determining the coefficients 
of each nodal point, the program solves the six orbital parameters 
for each arc. 

The orbit parameters were assumed subject to a priori 
constraints of 20.0 meters in position and .005 m/sec in velocity. 

The a priori standard error was assumed' to be 1.0 meter for all 
altimetrie measurements. 

During the initial preprocessing of the GEOS-3 data, 
certain characteristics were identified for editing criteria. The 
first editing level automatically examined altimetry measurements 
for gross errors that occurred from data handling procedures, such 
as tape parity errors and measurement identification problems. -The 
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second editing procedure was to examine the altimeter measurements for 
continuity and to eliminate abrupt point to point changes. The final 
editing is based on a three-sigma criteria when compared to a polynomial 
smoothing function. 

3.4.2 Results. The reduction of GEOS-3 altimetric measurements 
over the North Atlantic included data from 112 passes. The surface 
grid formed by these passes was dense in the western portion, and very 
sparse in the northern. It did not extend to the southeastern portion 
of the North Atlantic. Figure 22 illustrates the surface geometry of the 
GEOS-3 ground tracks. The figure also reflects the nodal selection 
used for this reduction. Due to the variations in data densities 
from one area to the other, a proper value for k {Equation Section 
3.2) had to be determined for representing the entire North Atlantic 
grid. It was found to be important that the choice of 5 = ka bear 
a balanced relationship to the typical spacing of the nodes. The 
more closely spaced the nodes, the smaller the logical choice for 6. 

This allows the surface in the vicinity of any given data point to 
be determined predominantly by those kernel functions of nearby nodes 
and yet prevents any one node from exerting total dominance. 

Figure 13 shows contours of the derived geoid undulations 
using the multiquadric function approach. In general, the entire 
geoid agrees favorably with the Marsh and Chang gravimetric geoid, 
especially in the western portion of the area where the data density 
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FIGURE 12. Calibration Area GEOS-3 Ground Tracks 
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FIGURE 13. Geoid Contour Map of the GEOS-3 Calibration Area 





is much greater. The area with less data, such as the northern 
part of this region, did not reflect the full geoid detail shown 
by the gravimetric geoid used for this comparison. The rms of 
the residuals of the altimetric measurements after the adjustment 
was 1-8 meters. 

Geoidal profiles derived from two passes of GEOS-3 
altimeter data have been compared with the Marsh detailed gravi- 
metric geoid. The geoidal comparisons are shown in Figures 14 
and 15. As can be noted in both figures, the relative shape 
agreement for both passes is excellent. 
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SECTION i} 


MATHEMATICAL SURFACE MODELS 

As stated earlier, the model originally used in SARRA 
for representing the oceanic geoid, was a derivation of the spheroidal 
multiquadric analysis developed by Hardy {1972). This model provides 
convenient flexibility when processing data collected over extremely 
small surface areas with a dense set of measurement data. It is 
independent of spherical harmonic representation of the geoid, thereby, 
eliminating the requirement to solve for high order and degree spherical 
harmonic coefficients. However, in the course of preliminary data 
processing of GEOS-3 altimetric data, it became evident that the 
observation data was very sparse in certain areas of the North Atlantic 
Ocean. Therefore, another model, namely the covariance function by 
Eeishmen and Moritz ms investigated and implemented into the SARRA 
computer program to provide more flexibility when processing sparse 
sets of data. 

4.1 The Covariance Function 

The detail derivation of the covariance function as applied 
in the SARRA computer program may be found in Appendix A of this re- 
port. The expression for the partial derivative matrix is defined 
here for the formation of normal equation. 
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The covariance function for geoid undulations is 
obtained by averaging the product N over the unit sphere, 

D(i|)) = (1) 

N and in this expression are the geoid undulations at any two 
points separated by the symbol M indicating th^e average over the 
unit sphere, namely 

M {•} = — (*)dcr, 


where a represents the surface of this sphere and da is the element 
of surface area. The covariance function may be expressed as {Eeiskcmen 


and Moritz i 196?)'. 

CO 



n=2 


n 


tcos 1 |;), 


( 2 ) 


Where 

V -2 -2 . 

n-o 

the quantities dfj are called "degree variances" for geoid undulation 

and A? are the corrections to the a priori coefficients of the 
nm — ■' — “ — 

spherical harmonic expansion of the geopotential. The parameter R 
is the mean radius of the earth and G is the mean gravity of the 
earth's surface. By ij;=0, we have from (2) by substitution 


D{0) 




n=2 


n 


( 2 ^) 
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The weighted sum of covariance functions between point 
and a number of selected nodal points Pj gives the undulation 
P^: 

where cy are the weights which must be determined and the 

spherical distance between and Py. 


= 


-I 


(4) 


4*2 The Spheriodal Multi quadric Model 

The expression for the geoidal model expressed in terms 
of the spheroidal multi quadric functions is the same equation (4) 
above. The only difference is the replacement of the covariance 
function D(i|j^y) with the "kernal function" of the form 




- Xyl + fea)2 + {|Y^ - Yyl + ka)^ + (|z - Z;l + 


U) 


Where 

cos (j)^ sin X^ 
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And 


Hr = Altimetry measurement of the geoid point 

r^ = Geocentric radius to the geoidal point Xg^, Yg^, Zg^ 

fea = An arbitrary fraction of the semi -major axis a 

X;, Y*, 1 - ~ arbitrarily specified nodes. 

J J J 

Xg^, Yg^, Zg^ = geocentric coordinates of the i.th 
sub-satellite point 
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SECTION 5 


FORMATION AND SOLUTION OF THE NORMAL EQUATIONS 

The observation equation for a single pass is taken from 
equation (13) of Appendix A and written as follows: 


V 


X. 





6(R/ - ril '] X 
6(X,...2)o J 


+ (R| - rf - ), (1) 


where 

= residual 

= initial estimate of the radial distance to the sub-satellite point 

d^ = small quantity to correct for non-geocentric direction of the 
altimetry measurement 

dc^ = coefficient matrix of nodes 

dx,...dz = coefficient matrix of orbit parameters 


dCi 

dx 
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Here the subscripts Z and j refer to the sub-satellite 
■point and the node of a single pass. Later the subscript k will 
be introduced to refer to the arc. The ‘‘kernal function", 
Dl’I'yyOs is selected by input to be either the covariance function 
or the multiquadric function to represent the geoid surface as 
described above (Section 4.1^ 4.2). 

Equation U) is written in terms of the matrix of partial 
derivatives for the altimetry, measurement point in the form. 


( 2 ) 


B; = 


_ 1 




Where 


Where 


= standard error of .the altimeter measurement 


B; taken from equation (1) as 


13 ) 


By = - 


.... 




n = number of nodal points 


Bx: = 3{R / - r_4_ 

3(X , . . . Z) Q 


( 4 ) 
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The parti als with respect to the orbital parameters have 
been refined by Bldka {1976} which introduced a faster convergence 
in the adjustment and in most cases removed the need for iterating 
the solution. The detailed description of the parameters in (4) 
may be found in Appendix A. 

The discrepancy term is defined by: 

= lie ~ n - + dz . (5) 

The normal equations for the arc are written as: 

% ( 6 ) 


where 


m 


Nfe = 


m 




6^ = correction vector 


and 


.th 


m = number of altimeter measurements in the k arc. 
The normal equations generated by this arc may be written 
in the matrix form 










) 

. 

J 


( 7 ) 


Introducing a priori orbital constraints for the arc 
leads to the following system of normal equations. 
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( 8 ) 






* 

■«fe % 

- 

k 


% 

4 + “fe 

V- 


• « 


_(‘4- •ci 5 _ 


Where 

llof^ = the inverse of the covariance matrix of the 
a priori values of the state vector. 

= the difference between the a priori values of 
the state vector and the values currently being 
employed as current approximations. Initially, 
i may be considered zero for the first iteration. 

The expansion of equation (8) is now introduced for the 
simultaneous reductions of all adjustable parameters for all arcs. 


’ N ■ 1 Ni Nj Nj 
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1 

Cc - *0)' V)s 


(9) 
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in which 


( 10 ) 


N = Z) "fe 

k=l 

3 

The solution of the normal equations ^9) involves n unknowns 
[one for each node) plus six for each arc. The number of arcs, s, 
in a solution may grow to several thousand, thereby, creating a large 
system with rank, r = n + 6s. The solution to this large system of 
normal equations is made practical by virtue of exploiting the patterned 
characteristics. Full details of derivation may be found in Brown [1958] 
or in Brown^ Trotter {1969). The computational steps follow. 

Compute the auxiliaries for the pass 



(N +' m3*\ 
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m,6) 

16 , 1 ) 

(6,6) 

(6.1) 
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As the and c.^ matrices are formed for each pass, they are summed 
into the matrices S and a by 


S 





and 


c. 


s 


L 


% 


(15) 


(16) 


The solution for the geoidal coefficients is computed by 
S-1 c 

(n,l)(ri,n)(n,l) 


(17) 


The solution 6^ for the adjusted corrections to the state vector 
for the arc is computed by 



(N 


I « • 


.-1 


(Cfe - Wfe ‘4 ) - Qfe Sfe- 


(18) 


• The error propagation for a given geoid point may be 
computed by 


a- . 

AyL 


= O: 


yCL 


■<V‘ 




(19) 


where the Bp partial derivative matrix is computed for the corres-- 
ponding 4>> ^ as in equation (3). 
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’ SECTION 6 


ANALYSIS OF RESIDUALS 


Careful pass-by pass examination of the plotted altimeter 
residuals can be made in order to uncover any localized systematic 
effects signifying insufficiently modelled geoid undulations. This 
is especially important to areas corresponding to known bathymetric 
and' surface features such as the Puerto Rico Trench, the Mid-Atlantic 
Ridge, sea mounts, etc. Where indicated from analysis of residuals, 
additional nodes can be introduced at appropriate locales. The 
process of nodal densi fication for more' detailed local definition 

of the geoid will exercise supplementary observations from any 
additional passes that may be available over that area. 

Altimeter residuals will reflect not only unmodelled 
geoid undulations, as just discussed, but also various quasi- 
systematic trends attributable to ocean dynamics such as tides, 
wind stress, swells, currents, etc. With the exception of currents, 
such phenomena are ephemeral or cyclic and thus, random over a 
sufficiently large number of passes. Except near shore lines, their 
amplitudes are generally less than one meter and will have only a 
slight influence on the geoid recovery. Although the residuals 
obtained following nodal densifi cation contain a wealth of infor- 
mation, it is beyond the scope of this contract to subject the 
residuals to such a detailed analysis. 
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This section will present the procedures that may be pur- 
sued in facilitating the ultimate residual analysis. As thousands 
of passes have been processed with the computer program SARRA over 
a course'of months or even years, a very stable data base of normal 
equations and coefficients of the geoid surface will be established. 

As new observation data is collected {es-pecially in regard to SEASAT A), 
one only has to form the normals of the new passes and sum them into 
the accumulated normals saved in the data base in order to update 
the solution. 

Once the geoid coefficients become sufficiently stable 
[meaning additional measurement makes insignificant changes to the 

I 

geoid} t we Bjay ‘treat the geoid coefficients as being perfectly 
known and process individual passes for detail residual analysis. 

The only adjustable parameters are the six state vector parameters 

♦ * * 

(jsfj Zj Zj J, z). This leads to a simple adjustment that could 
be processed on a mini -computer or a hardwired microprocessor 
requiring less than 1000 words of memory. It could even be envis- 
ioned as a real time monitoring processor with the appropriate com- 
munication and transmission interface. 

The solution of the single pass may be accomplished as 

follows; 

the observation equation is written as 
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(1) 


r = (1^ - '■'U R ” 

d (X,Y,Z,X,Y,Z) d Z 




where all parameters are identical to those described in Section 4. 

The matrix of partial derivations of the altimeter measurement 
is 

the discrepancy term 


^ d (X.Y,Z,X, Y,Z) 




The geoid parameter, r^, is computed from the geoid . 
coefficients, 6, equation (18, Section 5) saved. from. the data base {previous 
SAREA reduotion) and parti als, , of the geoid parameters for the 

measurement. The matrix is computed and formed by the same 
process as in Section 4.1, equation (2) or Section 4.2, equation (1)» 
depending on which model was used to represent the geoid surface in 
the master reduction. The solution for r^ is 


r^ = 6 (4) 

(1,1) (l,n) (n,l) 
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where 


n = number of nodes used in the master reduction 

and the parti als, are computed for all nodes in 
the master solution. 

The normal equations are written in matrix form as 


where 


and 


N 


6 


c 





B. 

-c 


c = 


6 = 


m 

I 

1=1 

d 
d 
d 
d 
d 
d 


•B-[ 


E. 

A, 


m = number of observations for the arc. 


(5) 


In order to make the above solution determinable for 6, the orbital 


constraints must be applied as in Section 5. 

Equation (5), above, is modified to accommodate orbital 
constraints in the form 


• • •• » » •• •■ •• 

(N+w)5 =c-w^ 


( 6 ) 


48 



Since the value for ^ is normally zero unless some correction is 
known in regard to the state vector, the solution for 6 becomes 


6 = (N + 0 ))-^ c 

( 6 , 1 .), ( 6 , 6 ). ( 6 , 1 ) 


(7) 


By substituting S into equation (1), we solve for the residuals, 
r^ for all points on the arc. These could be output to plot displays 
on a CRT type monitor. 

The above process requires minimal computer core requirements 
and the largest matrix required is (6,6) for the N and w matrix. The 
processing time would also be minimal. 
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SECTION 7 


FINAL REDUCTIONS OF THE NORTH ATLANTIC OCEANIC GEOID 

The major objective of this effort was to develop the 
appropriate computer software for reducing altimetric observations 
and to apply this software in the determination of the North Atlantic 
oceanic geoid. The initial coverage was expected to be similar 
to that depicted in Figure 1 and become more densified as more 
observations were made by GEOS-3. However, the availability of 
data did not provide the. desired coverage over the entire North 
Atlantic, nor the density in key areas outside the calibration 
area. This data set provides enough coverage {see Figia>e 16) for 
an extensive testing of the capabilities of the software system in 
detecting and modelling detail features. 

The final reductions are essentially extensions to the 
preliminary reductions presented in Section 3.4 

7.1 Unknowns and Constraints Used in the Reductions 

The same set of nodal points, data and orbital constraints 
were used for both the spheroidal multiquadric and the covariance 
function models. The number of nodes used was 47 and the number of 
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passes was Z36. This resulted to a total of 1463 unknowns {6 x 

passes + nodes) that were recovered simultaneously. The constraints 

exercised for the state vectors in terms of north, east and up were 

cr = 20. 0 metevs 

n 

CTg = 20. 0 meters 

= 20. 0 meters 

= .010 meters /second 
Og = .010 meters /second 

= .010 meters/second. 

These orbital constraints are chosen to be relaxed enough to accom- 
modate expected errors in the generated state vectors. The alti- 
metric observations were treated as subject to 1.0 meter error {one 
sigma) for the adjustments and error propagation. 

The selection of nodal positions are illustrated by Figure 16. 

7.2 Reduction Results 

The reduction of the GEOS-3 aitimetric measurements 
include data from 236 passes of GEOS-3 over the North Atlantic. 

The surface grid formed by the ground track of these passes was 
very dense in the western portion and sparse elsewhere. Figure 16 
illustrates the surface geometry of the grid. 
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FIGURE 16. Ground Track of Satellite Passes over the North Atlantic Ocean 





It has been experienced [vefevenoe 11) in previous appli- 
cations of the spheroidal multiquadric model that a proper choice 
of 6 = fea {EquaHon 1, Section 3.2) is important in keeping a bal- 
anced relationship to the typical spacing of the nodes. The constant 
6 controls the degree of correlation between nodes and as the value 
of 6 decreases, the correlation between nodes decreases (see refer- 
enoe l) . The present reductions varied in data density from one 
portion of the North Atlantic to another, thereby, making the choice 
for 5 a compromise between the dense arid sparse areas. This is prob- 
ably the primary weakness of the use of the multiquadric model. How- 
ever, provided uniform grid, adequate density and the proper selection 
of 6, the geoid surface can be measured to a high degree of local detail. 

Figures 17 and 18 represent the geoid contours from the re- 
ductions utilizing the covariance functions and multiquadric- models, 
respectively. The two contours generally agree to about one to two 
meters in the very dense portion of the North Atlantic. The less 
dense portion of the North Atlantic differs by as much as ten meters. 

When compared to the Marsh and Chang geometric geoid, it appears that 
the reduction using the covariance function agrees more favorably, 
although, the reduction using the multiquadric model showed better 
agreement in some of the more densely grided areas. 
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The geoid contours that are presented in Figures 17 and 18, 
are automatically generated by evaluating the nodal coefficients 
obtained from the adjustment for a set of coordinates destributed 
throughout the observed area. 

7.3 Residual and Profile Analysis 

Once the coefficients of the nodal points are determined, 
the resj duals are computed along each arc for evaluating the point 
to point fit of the observations to the adjustment. As stated in 
Section 6, these residuals reflect certain unmodelled surface fea- 
tures in addition to the random noise of the altimetric observations. 
No attempt will be made in this report to evaluate the meaning of 
these residuals since it is beyond' the scope of this contract. 

In addition to computing residuals, the SARRA computer 
programs computes geoid profiles along each arc. These profiles and 
the residuals are saved on a magnetic tape or disc during the SARRA 
adjustment and on option may be plotted for a graphic display. The 
residuals and geoid profiles of a few selected passes are presented 
in Figures 19 through 27 as examples. 
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FIGURE 19. Residuals from Pass 2S31413 
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FIGURE 20. 


Residuals from Pass 1222133 
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FIGURE 21. Residuals from Pass 3040601 







FIGURE 25. Geoid Profile along Pass 1222133 Using 
the Covariance Function Surface Model 
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FIGURE 26. Geoid Profile along Pass 3040601 

Using the Multiquadric Surface Model 
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SECTION 8 


PRE AND POST PROCESSING AUXILIARY PROGRAMS FOR 
SARRA REDUCTION OF SATELLITE ALTIMETRY DATA 

Several computer programs were developed under this contract 
to perform editing [ppe-pr*oeessoi’) prior to SARRA reduction and dis- 
play [plots) results after the SARRA reduction. The pre-processor 
served to edit raw altimetric measurements on both a gross error 
and on a statistical basis. In addition to editing the data, the 
pre-processor also provides a means of computing a mid-arc state 
vector, given the sub-satellite position. This was required 
to .process some data provided by NASA where only tj>, X and h were 
given. The plot program serves 

(a) to plot satellite ground tracks of 
available data^ 

^b) plot contour map of resulting geoid, 

[c) geoid profiles along a given are, 
and 

(d) residuals from the adjustment. 

8.1 Pre-processor and Data Noise Filter 

Pre-processing and editing [PRBP] of GEOS-3 altimetric 
data is designed as a three-level automatic computer editing effort. 
The raw data is fiVst edited for "gross" errors. The second step 
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is to examine the data for continuity and eliminate abrupt changes. 

The third level of editing rejects data, either smoothed or unsmoothed, 
that have standard deviations greater than some specified criterion. 

The editing is accomplished in three steps, 

(l) ootnpap-ing aZtimetrio measuvements with a 

veasonabte eospeoted range of 'possi-bte heights 
(k) for gross error such as parity errors; 

(.2) continuity is detected by comparing the change 

in the measured altimetry (Ah) bel^een successive 
data points with the corresponding change in time 
(At), This is a measure of the altimetry rate 
(h) which is a smooth function. If the altimetry 
rate exceeds the maximum value y the point is 
refected. The maximum for h may be estimated 
emperically from several sets of GEOS-C altimetry 
data. 

(3} The third level of editing refects data based 
upon an input sigma criterion. The criterion 
will be based upon realistic expected accuracy 
of the GEOS-3 altimeter measurements. 

The criterion chosen in step 3 should be, a little relaxed 
to make sure that points are not rejected that truly reflect some 
surface details that are normal sea state variations. 

In addition to editing measurement data, the pre-processor 
can compute a mid-arc state vector when the ephemeris is provided 
with the alti metric data tape in the form of latitude, longitude 
and height. The following computational steps are implemented in 
the pre-processor. 

The procedure followed in computing a mid-arc state 
vector for each pass, was to compute Cartesian coordinates from 
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the provided geographic coordinates. These coordinates were 
computed by 


N = a/(l - sin^ <^) 


1 

2 


"x" 


(N+h) cos cos A 

Y 

= 

(N+h) cos 0 sin A 

.z. 


_ [N{T-g 2) .+ h] sin 


Where a = the earth semi-major axis 
e = eccentricity of earth 


A fifth order polynomial fit to these coordinates was performed 
in order to reduce the effect of the above truncation error and 
to provide a means for computing velocity components. 

The general polynomial expression is given by 


The equation written in matrix form for all i is expressed by: 
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Substituting C for the X. vector, B for the t’J matrix and S for 

A- 

the a- vector, yields 
C = BS . 

The least squares solution {with vnit weight matrix) for the 
coefficient vector is 

S = (b' B)^ b' C . 

til 

The n order polynomial equation is evaluated at mid^arc time, 
where 

fo ■ - t^) / 2 , 

The velocity component at X, is computed from the first 
derivative of position by 

X = ^ = =1 + 2 Vo -"3 ajiJ + na i" . 

This procedure is repeated for the Y, Y and Z, Z 
components. All altimeter observation times (t^) are initial- 
ized with respect to mid-arc by 

t . = t. ~ . 

AX) A, O' 

The velocity components corrected for the earth rotation are 
computed from 

*0 = *0 - ’’o 

y„ = Yg + TXg 
• 

to obtain inertial components (¥ = earth rotation rate). 
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The resulting state vector representing the orbit at 
mid-arc {Xq, Zq, yo and z^) is input to the SARRA program 
with the appropriate constraints. 


8.2 Post SARRA Plot Programs 


The plot program provides the capability of producing satellite 
ground track plots, geoid contours and either geoid profiles or 
residual profiles. Additionally, the land features superimposes on 
the ground track and the contours.’ Initially, the land outline was 
read from maps and input into the program plots. This was later 
replaced at AFGL by a magnetic tape containing the digitized shore- 
line coordinates. Each plot is described as follows: 

(a) Ground Tfuidz Vtot . This program plots the 
ground traek of satellite passes over' the 
North Atlantic from a given sub-satellite 
latitude and longtitude. Additional paral- 
lels and meridians may he selected with 
program input parameters. - The ground track 
plots may be contained yyithin the North 
Atlantic data houndaries as shown in 
Figure 28. 

(b) CovitouiA. VZott. Ooivif-GIA^S Oj^ G'l' ther The 
standard deviations or 'The actual geoid 
heights may he plotted. The standard 
deviations are computed from the covari- 
ance matrix from thie SABRA reductions. 

The geoid heights are computed from The 
surface coefficients saved from the 
SARRA solution.. In either case. The 
computations ore performed for a grid 
of surface points. -The generated 

grid is contained within the hound- 
aries of the observed geoid and the 
houndaries of the North Atlantic. 
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(c) R&6Zduat6 . 2%e dltimetev observation 

residuals (me.a&uA.&d mi.yuL6 c.omput 2 .d 
hzight) obtained from the SARRA 
reduetions are plotted for the pur- 
pose of visually reviewing measure- 
ment oharacteristios . The residual 
, profile plot provides immediate' dis- 
play of altimeter measurements noise 
and such oceanic details as sea state 
variations . 


71 




APPENDIX A 


The use of the Covariance Function In the Short Arc Reduction 
of Radar Altimetry Program 

The following derivation of the covariance function and its 
application to the modelling of geoid surface in the SARRA computer pro- 
gram was developed by Dr. Georges Blaha of DBA Systems, Inc. under a 
special investigation. The results of Blaha's investigation is re- 
produced as an appendix for detailed reference material regarding 
Section 4 of this report. 

Covariance Function for Geoid Undulations 

The covariance function for geoid undulations (N) is derived 

from the basic formula 
» n 

, n *=2 m=o 

where 

R - is the mean earth radius (R sj 6371 km), 

Y - is the mean value of gravity (y - 979.8 gal), 

"I", X - are the geocentric latitude and longitude of 
the point associated with N: 

the other overbars indicate that we are dealing with "fully normalized" 
harmonics. We have (m / 0): 

\o \ f) = 42 n + 1 (sin 

5) = -^2 C^n i). 
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Pjj (sin qp) being called Legendre's polynomials and 

associated Legendre functions, respectively. Furthermore, 




the coefficients C , C , S are sometimes called "conventional C*s 

no nnr nm 

and S’s", 


It has to be emphasized that et^uation (1) is valid for the reference 
ellipsoid having the same mass and the same potential as the geoid. This 
condition is fulfilled by the "mean earth ellipsoid" which in theory 
shares with the actual earth two additional parameters (iu = rotation rate 
and J 2 = -C 20 )• The pertinent constants of the mean earth ellipsoid in 
the Geodetic Reference System (GRS) 1967 are: 


c* 

so 

= - 1.0827 X 10" , 

CAa) 

C40 

= + 2. 37 X 10 ",^ 

(4b) 

Cto 

- 6. 1 X 10 ; 

(4c) 
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the other C* are essentially zero (already C* is very small and is 
no 60 

sometimes neglected). 

In what follows we shall be dealing with overbarred coefficients 
related to their conventional counterparts by the factors shown in the 


equations (3) . Adopting the mean earth ellipsoid as the reference 


ellipsoid^ we thus have in agreement with the equations (4) : 


C*o = -.484194 X 10 “3, 

Cjo = +.790 X 10“^ 

^ = -.17 X 10"®. 

60 

In analogy with previous notations (conventional case), we write 


A C 20 

— C 20 “ 


= ■^0 - 

A C 


60 

60 



(5a) 

(5b) 

(5c) 

(6a) 

(6b) 

(6c) 


For all the other A we have essentially (even if m = 0) : 

AC = C ; 
nm nm 


(7a) 


all A S are in fact S themselves 
nm nm 

AS = S . 
nm nm 


namely 


(7b) 


The covariance function D (ijr) for geoid undulations is obtained by 
averaging the product N N‘ over the unit sphere, 

D (i|f) = M (n N‘3 J (8) 

N and N’ in this expression are the geoid undulations at any two points 
separated by the spherical distance i|r and the symbol M indicates the 
average over the unit sphere, namely 
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I! 

a 


where O represents the surface of this sphere and dQ is the element of 
surface area. It can be shown that the covariance function in (8) may 
be expressed from (1) as follows: 


D (t) 




^n ^n (*^0® 


(9) 


where 



(n- 1) = 


n 

I 


(A + A .S^ ) ; 
nm nm ’ 


m 


( 10 ) 


the quantities are called "degree variances" 
If ^ = 0, we have from (9) : 


D(0> 




n — 2 


d . 

n 


for geoid undulations. 


(9’) 


In the program SARRA, the covariance function serves to express 
the geoid undulation (%) at a geoidal point as a function of 
parameters Cj at selected nodes P^: 

Z * 

D (llTij) Cj , (11) 

J 



adjustment model in SARRA can be written as 


The overall 


Hi = Ri - Tf + di. 
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where 


- is the radial distance (from the coordinate origin 0 ) to 
the satellite point 

r^ - is the radial distance to the subsatellite point 

H. - is the satellite altimetrv observation (distance S. P.). 
1 '11 


The small unadjustable quantity accounts for the fact that 0, and 
Sj do not lie on a straight line; it is not needed explicitly. The 
parameters c^ from ( 11 ) enter into the adjustment process through the 
model equation 


r j - + Nj, 


where r^ is the radial distance to the reference ellipsoid (it is computed 
f rom (p^) . The remaining parameters are the state vector components; 
these parameters are always weighted. The adjustment model in SARRA may 
thus be written symbolically as follows : 

(state vector parameters) - (c^ parameters) - r ’ ( 9 ) + d . , ( 12 ) 

When performing partial differentiation of this equation with respect to 
the parameters^ it should be kept in mind that cp in (cp) depends to a 
certain extent on the state vector parameters, (See reference 10 for a 
detailed discussion.) 


g2 _ t 

“ a/[l sin^ j s and being the usual ellipsoidal parameters 
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xn oroer xuna xiie e4,u^uxuu£3^ lluul we ueauce 


' S‘(C.. par. ) 

J o 


^ (R. - r!) 

r' + d + 'r = 

i i o (state vector par.) 

o 

d (Cj par.). 


d (state vec. par.) 


where is the residual, is the observed altimetry and "o" denotes 

the approximate (initial) values of parameters or functions of such 

« • • 

parameters. The state vector parameters are X, Y, Z, X, Y, Z in the 
"Earth Fixed" (E.F.) coordinate system; the adjustment in SARRA is 
effectuated in this system. From (11) we have 


S N, 

1 

d (c. par.) 

J o 


— {D h • • • ] 5 


due to the linearity in (11), one can take 


c” = 0 


and thus 


N? = 0, 


which is in fact used in SARRA. 


The observation eauation is then written as follows; 


Vf = ^ - [D (tj_^),..., D (tj.j),...l 


dci 


S (X,..,Z)^ 


X 




dX 

di 


+ (R° - r| - hJ + d^), 


(13) 
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With self-explanatory notations this equation is presented as 


Vi = ^i y + 1^, 


(14a) 


where 




= R° - + d . . 

i i 1 r 


(I4h) 

(14c) 


The term 1^ is sometimes called "discrepancy term". The quantity r. 
can either be computed as in (14c) ^ or it can be computed as the distance 
OP^ from the coordinates of the point depicted in Figure 1, This is , 



By joining together individual observation equations along one short 
arc, the resulting set can be written in matrix notations as follows: 


V = A )f + L. 

Similar procedure would apply for other arcs in an adjustment; each vector 
X would now contain a different set of state vector parameters (besides 
the common set of Cj parameters) . All these sets would finally be joined 
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D(i') 

2 

(METERS ) 



FIGURE 2: COVARIANCE FUNCTION FOR GEOID UNDULATIONS 

ASSOCIATED WITH DEGREES TRUNCATION N=2, N=6, 
AND N=21 
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into one large system of observation equations, and adjusted. However, 
the main purpose of this brief exposition of the SAREA adjustment model 
has been to show the role of the covariance function D (i|/) . We have seen 
how to compute this function from potential coefficients. We shall next 
discuss some practical problems associated with its application. 

If a set of potential coefficients (C's and "S's) is given complete 
through a certain degree n = N, the covariance function D (\|r) may be 
constructed for varying values )jr (between 0 and 2 tt ) according to the 
equations (9')^ (9). In fact. Figure 2 depicts three cases of the 
covariance function D (ifr) ; they correspond to the truncations N = 2, 

N = 6, and N = 21, associated with a reasonable set of coefficients 
complete through the degree and order (21,21). The figure reveals that 
the greatest contribution to the value of the covariance function comes 
from the lower degree and order coefficients; there is a small difference 
between D (i}r) for N = 6 and N = 21 and there would be practically no 
difference if_N = 15 and N — 21 were considered. On the other hand, 
it appears that if, for example, a given (6,6) model were to replace 
the reference ellipsoid, the corresponding covariance function of new 
’’undulations** would exhibit sharp and distinct features which would vary 
substantially for different degrees of truncation (up to a relatively 
n = N) , It is felt that such an approach would be much more 
sensitive to local geoidal features than is the use of the covariance 
function D (ijf) . In contrast to using D (i|f), such a ”modified covariance 
function” would greatly reduce the influence of distant nodal points (e.g. 
%fhen ^ >30®) on the value of the "undulation" at observation points. 
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